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ABSTRACT 

The results of a spectral analysis, using XMM-Newton and Chandra data of the brightest ultra luminous X- 
ray source in the nearby galaxy M82, are presented. The spectrum of M82 X-l, was found to be unusually 
hard ( photon spectral index T w 1 ) with a sharp cutoff at w 6 keV. Disk black body emission model requires 
a nonphysically high temperature. Instead, the spectrum is better described, with a lower reduced % 2 , as 
emission due to nearly saturated Comptonization of photons in an optically thick (r w 10-30, depending on 
the geometry) plasma having a temperature kT « 2 keV. This is in contrast to the high energy spectra of other 
black hole systems, which are relatively steeper (T > 1 .5) and hence are modeled as un-saturated thermal and/or 
non-thermal Comptonization of soft photons, in an optically thin (r ps 1) high temperature plasma. An Iron 
line emission which is marginally resolved (er ~ 0.2 keV) is required to fit the data. We argue that the standard 
geometries for the X-ray producing region, which are an optically thin inner disk or an uniform/patchy corona 
on top of a cold disk, are not applicable to this source. Alternatively, the geometry of the X-ray producing region 
could be a large sphere surrounding a cold accretion disk or an optically thick inner disk region which cools by 
bremsstrahlung self-Comptonization. For the latter scenario, such an inner disk region, whose effective optical 
depth to absorption is less than unity, is expected in the standard accretion disk theory for near Eddington 
accretion rates. 

Subject headings: Galaxies: general -galaxies: individual (M82)- X-rays: binaries 



1. INTRODUCTION 

Chandra observations of nearby galaxies have detected 
several non-nuclear X-ray point sources, some of which 
have isotropic luminosities > 10 39 ergs/s and are called Ul- 
tra Luminous X-ray sou rces, ULX ( Matsumoto et al. 1 I200T1 
Zezas & Fabbiano 2002). These results have confirmed and 
extended the earlier evi dence for the existence of such sources 
by the Einstein survey (Fabbiano 1989; Fabb iano & Trincheril 
119871) an d the ROSAT detectio n of ~ 100 ULX in nearby 
galaxies llColbert & Ptakl 120021) . Systematic studies of the 
Chandra observations have revealed that ULX exist in galax- 
ies of different morphological types. Since, by definition, 
ULX emit radiation at a rate larger than the Eddington lu- 
minosity for a ten-solar mass black hole, they are believed 
to harbor a black hole of mass 5OM <M< 1O 4 M , where 
the upper limit is constrained by the fact that a more mas- 
sive black hole would have settled into the nucleus due to 
dynamical friction (Kaar et et al. 1 12001). If this is true, then 
the differences and similarities between the spectral proper- 
ties of such intermediate mass black hole systems (IMBH) 
with those of Galactic black holes (GBH) and Active Galac- 
tic Nuclei (AGN), are expected to enhance our understand- 
ing of black hole systems in general. Alternatively, the ULX 
sources may be systems emitting at super-Eddington lumi- 
nosities (Begelman 2002), or their emission may be beamed 
from a geometrica lly thick accretion disk ( [Reynolds et al. I 
119971 iKingl 120021) . It has been argued that in the latter 
case, such thick "funnel" shaped disks can enh ance the ob- 
served flux by just a factor of few (Misra & Sriram 2003). In 
any case, the com plexities involved in explaining their origin 
JPortegies. Zwart & M cMillianl l2002t iTaniguchi efaTlteOOrt 
Madau & Reesl 120011) and the inferred high accretion rates 
( King et al. 1 1200 llT make these enigmatic sources an impor- 
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tant class of X-ray systems, whose dominant radiative mech- 
anism and geometry, needs to be studied and constrained. 

The spectra of ULX are expected to be different from GBH 
and AGN, if the black hole is an IMBH and/or if their lu- 
minosity exceeds the Eddington limit. Earlier analysis of 
ASCA data revealed that the spectra of some ULX can be 
described as radiation from a black-body emitti ng disk with 
an ad ditional high energy power-law (Colbert & Mushotzkv 
[19991 which is also the standard radiative model for the 
high-soft spectral state of Galactic black hole systems. The 
inferred inner disk temperature from these initial analysis, 
was found to be too hig h (kT,„ ~ 1.1 — 1.8 keV) for an ac- 
cretin g IMBH (IColbert & Mushotz k\ll999t rMakishim a et alJ 
2000). However, recent analysis of XMM-Newton data have 
revealed the presence of a cool accretion dis k component 
Qdjn ~ 0.1-0.5 keV) in NGC 1313 X-l, X-2 (iMiller et al. I 

2003) and M81 X-9 ( IMiller. Fabian & Milled 12004 which 
reaffirmed that these systems may indeed harbor an IMBH. 
The spectra of some other ULX can be represented by a sin- 
gle power-law suggesting a similarity with the low-hard state 
of GBH, where the disk emission is either absent or weak. 
In fact, transitions between the two spectral states, similar 
to those seen in GB H, have been repor ted in the ULX NGC 
1313 X- 1 (IColbert & Mushotzkvl 1 1 9991) and two sources in 
IC342 ( Kubota et al. Il200l | ). Analysis of Chandra data for a 
large number of ULX (e.g. IColbert et al. l2004USwartz et al. I 

2004) have revealed that the spectra of most of these sources 
can be represented by a power-law with a typical photon in- 
dex, r pa 1.6-2.0. This spectral index is si milar to those of 
GBH i n the hard state and Seyfert 1 AGN (Zdziarski et al. 
1996a). However, most of these ULX are too faint for any 
further detailed spectral analysis. 

One of the most intriguing ULX is the brightest X-ray 
source in the starburst galaxy M82 that is located ~ 170pc 
away from the dynamical center. The isotropic luminosity 
of this source, M82 X-l, can exceed 10 41 ergs s -1 indicat- 
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ing that either it is radiating well above the Eddington limit, 
or it harbors a black hole > lOOMo, even if some moderate 
beaming (a factor of few) is taken into account. ASCA ob- 
servation have shown that this source is highly variable with 
its luminosity in 0.5 - 1 keV band changing from 4.5 x 10 40 
to 1.6 x 10 41 ergs s _1 jPtak and Griffithsl fl999). The pres- 
ence of a 54 mHz QPO in its X-ray light curve suggests 
that the source is a compa ct object and not highly beamed 
JStrohmaver & M ushotzkvl l2003l) . If this QPO is related to 
the 0.8 - 3 Hz QPO observed in galactic sources, it would im- 
ply that its black hole mass, M B h ~ 100-300M©. Its high 
luminosity and relative proximity makes this source an ideal 
ULX whose spectra c an be analyzed in detail. 

Swartz et al. (2004) have analyzed and fitted a large num- 
ber of ULX observed by Chandra including M82 X-l. Af- 
ter correcting for pileup, they report that the 1999 Chan- 
dra observation of M82 X-l can be fitted by a power-law 
with sp ectral index T = l-17+g;|| and N H = 0.98^ x 10 22 
cm -2 . Strohmaver & Mushotzkv (2003) mention that the 
XMM-Newton data can be fitted both by a disk black body 
emission with an unphysically high temperature of ~ 4 keV 
and a Comptonization model, but they do not quantify the 
parameter values obtained nor discus the physical in terpre- 
tation of the fits. Recently. iFiorito & Titarchukl (12004). fitted 
the same XMM-Newton data with the bulk Comptonization 
model and obtained a power-law index r ~ 2 indicating that 
this source is similar to the soft/high state of black hole bi- 
naries. However, as shown in this work, the column density 
required for such a fit, Nh > 5 x 10 22 cm -2 , is not consis- 
tent with some of the Chandra observations. Note, that the 
Chandra ACIS data essentially measures the spectrum in the 
0.5-5 keV, since the effective area of the detector decreases 
rapidly for energies greater than 5 keV. On the other hand, the 
XMM-Newton EPIC data measures the spectrum in the 3-8 
keV range, since below 3 keV the counts have contribution 
from the diffuse background. This complementary nature of 
the two types of observations makes it imperative to take into 
account both of them in any spectral analysis. In particular 
the Chandra observations give a better measure of the neutral 
column density, while the detailed higher energy shape of the 
spectrum can be constrained by XMM-Newton observations. 
Since the source is known to be variable, in principle, any 
spectral analysis should be based only on simultaneous obser- 
vations, which currently do not exist. However, here we take 
a pragmatic approach and use the non-simultaneous Chan- 
dra observations of the source, to test whether the parameters 
obtained are consistent with the results of the XMM-Newton 
analysis. Our motivation is to obtain the simplest physical ra- 
diative model for the source, which in turn, will highlight any 
differences between this ULX and other black hole systems. 

2. OBSERVATIONS AND DATA REDUCTION 

M82 was observed five times between September 1999 and 
June 2002 by Chandra. For three of these observations the 
pileup fraction for M82 X-l, is > 70% and hence no reliable 
analysis can be undertaken. For the other two ( whose ob- 
servation ID numbers and dates are tabulated in Table 1) the 
pileup fraction is around 25%, and hence only these two were 
analyzed. ACIS spectra of the source were extracted from 
circular regions of ~ 2" radius. Background spectra were ex- 
tracted from nearby source free regions. The analysis was 
done using the software, CIAO v3.2.1 and the latest calibra- 
tion files available till October 2005 were used to create the 



TABLE 1 
Chandra Data Analysis 



ObsID 


Date 


r 


N H 


C/s 


L 39 


PF (%) 


1302 
2933 


1999-09-20 
2002-06-18 


-Q-2 
11 -0.3 


, 2 +o.2 

1 1^ 

i - i -0.1 


0.12 
0.13 


4.8 

5.0 


25 
26 



NOTE. — Chandra data observational ID, date of observation, best fit 
power-law index (T), best fit column density in 10 22 cnT 2 (Nh), Count rate 
in the energy range 0.5-7.0 keV (C/s), estimated absorbed luminosity (using a 
distance D = 3.9 Mpc) for the same energy range in unit of 10 3 ' ergs/s (L39) 
and pileup fraction in percentage (PF) 



spectral response. 

XMM-Newton observed M82 for a total elapsed time of 
~ 30 ks on 2001 May 5. Both the EPIC instruments PN 
and MOS were operated in Prime-Full- Window mode using 
the medium filter. The PN events were processed using the 
XMM-Newton analysis software (SAS v6.0) and calibration 
data base available in August 2004. PN events corresponding 
to patterns 0-4 (single and double) were used for the spec- 
tral analysis. The PN spectrum was extracted using a circular 
region of 10" radius centered at the source position and the 
background spectrum was extracted from a source free (based 
on the Chandra image) circular region. The events were 
checked for the pile-up effect which was found not to be sig- 
nificant. There are two Chandra detected X-ray point sources 
located less than 13" from M82 X-l. To check whether these 
sources effect the analysis, we also extracted the PN spectrum 
from a smaller 5" region and performed the same spectral fit- 
ting as described in the next section, although the response 
for such a small extraction radius may not be accurate. There 
were no qualitative difference in the results of the spectral 
analysis. 

3. SPECTRAL ANALYSIS 

Since the Chandra observations are affected by pileup, the 
effect was taken into account by using the "jdpileup" model 
of the Chandra spectral fitting software Sherpa. We fitted 
Chandra data in the energy range 0.5-7.0 keV with "jdpileup" 
model convolved with an absorbed powerlaw. The best fit col- 
umn density Nh, the photon spectral index T, pileup fraction, 
count rate and estimated luminosity are tabulated in Table 1 . 
Since the pileup fraction is significant, fitting with a more de- 
tailed model than an absorbed power-law is not warranted. 
The qualitative inference that can be obtained from this analy- 
sis is that column density Nh ~ 10 22 cm -2 and that the spectra 
have a steep spectral index, T < 2. 

The XMM-Newton spectrum below 3 keV is expected to be 
affected by the diffuse X-ray background, hence only data in 
the energy range 3.0- 10.0 keV was used for spectral fitting. 
Emission f rom a black body em itting disk, the XSPEC model 
"diskbb" ( Mits uda et al. 11 984). does provide a formal fit (Ta- 
ble 2), but the required maximum temperature of the disk 
Ti n ^3.3 keV is unphysical and would require a large color 
factor correction ~ 10 for a disk radiating at the Eddington 
limit around a 100M Q black hole. Hence, a more plausible ex- 
planation is that the spectrum is due to Comptoniz ation. We fi t 
the data using the XSPEC model "Comptt" (Tita rchukl 1994 . 
which describes Comptonization of soft photons by a ther- 
mal distribution of electrons. There are two different sets of 
spectral parameters which can describe the data. The first set 
(second column of Table 2), has a high temperature, kT w 35 
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TABLE 2 

Best fit model parameters for XMM data 





U1SKDD 


fnmnTT 


TnmnTT 


Thcomp 


Bcomp 


nH 


O Q7+0.51 


c ca+0.46 
J ""-0.78 


2 47+ 2 - 12 


2 62+ L42 
Z - Dz -1.40 


9 77+1.40 
z - ' '-2.35 


m, 


-i -1O+0.13 
J - Jz -0.10 










kT 




3 =+30 


2 24 +0.28 
-0.22 


2 ^4+0.25 
z -> -0.12 


2 46+ ' 58 
ZTO -0.34 






0.30+jj;^ 


- n4 n 37 

- 3 °i8:ls 


0.l3*;f| 




T 




1 lfi+0 58 
llo -0.30 


11.10i|18 






T LE 








i 41+0.15 

1 -^ 1 -0.16 


i 4Q+0.23 
1 -^'-0.20 


Ei 


6 64+ ' 08 


6.4? 


6 67+ 08 
aD '-0.13 


6 67+ ' 07 
°- D/ -0.07 


6 67+0 .08 
0.0 /_ 13 


<*l 


u - z/ -0.08 




19+ 011 


io+o-'s 

uly -0.08 


02 ,+0.15 
uzl -0.09 


3a] 


172/176 


178/175 


171/174 


171/174 


171/175 


NOTE. — Diskbb: Disk black body emission model. CompTT: Comp- 
tonization model in Xspec. Thcomp: Comptonization model using the 
low energy spectral index Tle instead of r as a parameter. Bcomp: 



Bremsstrahlung Self-Comptonization model. For all models a Gaussian line 
at E\ and width ct; has been included. 
f E, fixed at 6.4 keV. 



keV, and low optical depth r«l, Such a spectrum, which 
is due to unsaturated Comptonization of external photons, is 
similar to the ones observed in other black hole s ystems and to 
the spectral fit obtained by Fiorito & Titarchuk (2004), using 
the bulk Comptonization model for the same data set. How- 
ever, the required column density nH w 5.5 x 10 22 cm -2 is 
larger than the that estimated from Chandra observations (Ta- 
ble 1). The second set of parameters (third column of Table 
2), has a low temperature, kT w 2 keV, and high optical depth 
t k, 10. This spectra, which is due to saturated Comptoniza- 
tion, gives a lower % 2 as compared to the unsaturated case and 
the required column density (which has large error estimate) 
is similar to the Chandra value. An Iron line emission, which 
is required for both the fits, turns out to be a narrow one at 
the expected centroid energy of w 6.7 keV (for ionized Iron) 
for the saturated spectrum fit. On the other hand, for the un- 
saturated spectrum, the required line energy is unphysically 
small at as 4.5 keV and has a large width of w 2 keV. Fixing 
the centroid energy to 6.4 keV (as reported in Table 1) still re- 
quires a large width of sa 1 keV. Thus, taking into account the 
smaller x 2 value, compatible column density with Chandra 
observations and the simpler Iron line profile, we consider the 
saturated Comptonization spectra as a better fit to the data. 

We confirm this result by also fitti ng the data to t he thermal 
Comptonization model described by ( Zdziars ki et al. 11 996b). 
named "Thcomp", which takes into account relativistic cor- 
rections (Table 2). Instead of the optical depth r, "Thcomp" 
uses the geometry independent low energy spectral index Tle 
as a parameter. The optical depth can then be derived from 
Tie and electron temperature, giving t sa 12 (sa 20) for a slab 
(spherical) geometry. The spectrum can also be interpreted 
as emission from a thermal plasma which Comptonizies in- 
ternally produced bremsstrahlung photons, instead of external 
low energy ones. Such a model, named "Bcomp", fits the data 
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FIG. 1 . — EFe unfolded spectrum of M82 X-l, for the soft photon Comp- 
tonization model. The filled circles are XMM-Newton data points which have 
been rebinned for clarity. The solid line is the unabsorbed model (while 
dashed line is the absorbed model) spectrum corresponding to the best fit 
parameters given in Table 2. 



equally well. Thus the fitting is not sensitive to the spectral 
shape of the input seeds photons used in the model. 

The Compton Y-parameter (y = AkTr 2 /m e c 2 ) is a mea- 
sure of the degree of Comptonization that an average photon 
undergoes. From the spectral fit parameters, y sa 3, which 
implies that the Comptonization process is nearly saturated. 
This is evident in Figure 1, which shows the unfolded EFe 
spectrum for the Comptonization model. The peak of the 
spectrum is at sa 6 keV sa 3kT, as expected for saturated 
Comptonization. This is in contrast to un-saturated Comp- 
tonization spectra which peak a t w kT. As reported by 
Strohmaver& Mushotzkv (2003) an Iron line feature is de- 
tectable in the spectrum. 

4. DISCUSSION 

The spectrum of M82 X-l is unusually hard and hence 
the simplest radiative model which can describe it is satu- 
rated thermal Comptonization. This is in contrast, to other 
black hole systems, where the spectra is generally mod- 
eled as un-saturated thermal/non-thermal Comptonization. 
The spect ra of other ULX are steeper with photon index , 
T ~ 1.7 tMrneret al. 1 12003 iMiller. Fabian & Milled l200l 
iMiller et al. 11200 4) similar to those of AGN and Galactic 
black hole systems. Thus either, M82 X-l is different from 
these systems or at least it was in a rare and unique spectral 
state during the XMM-Newton observation. 

The spectral fit parameters obtained imply that the geom- 
etry of the X-ray producing region for M82 X-l, is not the 
same as that of other black hole systems. The large optical 
depth t >10 would prevent the entry of any external soft pho- 
tons into the thermal plasma. Hence, the standard geome- 
try of an hot inne r disk, which Co mptonizi es photon s from 
an outer cold disk ( Shapiro, Liahtman & Eardlevl 19771) is not 
applicable to this source. Similarly, soft photons will not be 
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able t o enter an optically thic k patchy corona or active re- 
gion JPoutanen & Fabian 1999) on top of a cold disk. On the 
other hand, the soft photons from a cold accretion disk would 
diffuse out thro ugh an optically thick corona if it uniformly 
covers the disk (Liang & Prica ll977ft . However, such mod- 
els are applicable only when the Compton amplification fac- 
tor, A ~ 1 . Otherwise the Comptonized photons will impinge 
back on the disk and he at the disk to the corona temperature 
dZdziarski et al. Il996bl) . Since for saturated Comptonization 
the amplification factor A sa T /T s w 10, the uniform corona 
model is also ruled out. A possible geometry is that the in- 
ner region of an accretion disk is surrounded uniformly by a 
thermal plasma sphere. If the sphere is significantly larger 
than the soft photon producing disk, only a small fraction of 
the Comptonized photons will impinge back on the disk and 
energy balance can then be maintained. Interestingly, the fit- 
ted parameters suggest that the temperature of the accretion 
disk, which can be identified as the soft photon source tem- 
perature, T s ~ 0.3 keV, is expected for a disk around a large 
1OOM black hole. However, it is not clear whether such a 
large sphere would be dynamically stable or even realizable 
in the presence of a black hole. 

Another radiative model for the source is that the thermal 
plasma Comptonizies photons that are internally produced 
by bremsstrahlung. In th e standard accretion disk model 
(Shakura & Sunvaev 1973), the disk is assumed to be effec- 
tively optical thick to absorption and hence the disk radiates 
as a black body locally. However, at high accretion rates, the 
inner region of the disk may become optically thin to absorp- 
tion and the standard disk theory is no longer consistent. This 
transition occurs roughly when r* = y/T es T a b s s» 1, which cor- 
responds to a transition radius 



r tr w 40 r s a 



34/93 



M 



IOOMt 



2/93 



M 



M, 



Edd 



64/93 



(1) 



where r s = 2GM/c 2 is the Schwarzschild radius, M is the 
mass of the black hole, M is the accretion rate, M^d = 
LEdd/0-lc 2 , is the accretion rate corresponding to Edding- 
ton luminosity and a is the dimensionless viscosity parameter. 



For radii < r tr , the disk cools inefficiently by bremsstrahlung 
Self Comptonization (BSC) and hence is hotter than a disk 
described by the standard theory. In this region, the tem- 
perature increases rapidly with decreasing radii. In certain 
cases, especially when the effectively optically thin region is 
an extended one, the temperature can vary from «0.1 keV 
at the transition radius to nearly ~ 50 keV in the inner most 
radii. The resultant spectrum, which is a sum of local BSC 
spectra, can have a power-law form and such a scenario has 
been put forward as an alternate model to the standard un- 
saturated Compt onization interpretation of t he high energy 
spectra of AGN JMaraschi & Molendill200il) and X-ray bi- 
naries (Misra. Chitnis & Melial ll998l) . Since the temperature 
increases inwards, the high energy photons are preferentially 
produced in the inner region of the disk as compared to the 
lower energy ones. Thus the observed time-lag between dif- 
ferent energy bands (Nowak et al. 1999) can be explained as 
waves traveling from the outer to the in ner regions of the disk 
( Misra 200(l lUttlev & McHardvfcOOl . For M82 X- 1 , the sat- 
urated spectrum could arise from such an effectively optically 
thin region. However, here the temperature variation is more 
modest, with a maximum of ~ 2 keV, indicating that the re- 
gion is not an extended one. Detailed spectral prediction from 
such a model depends on the actual temperature profile, and 
the present data is not adequate to unambiguously constrain 
the physical characteristics of the region. Nevertheless, it is 
interesting to note that M82 X-l is probably emitting close to 
its Eddington limit, even if it is a 1OOM black hole, and at 
such high accretion rates, the inner region of the standard disk 
is expected to become effectively optically thin. 

The results obtained in this work need to be confirmed by 
simultaneous Chandra and XMM-Newton observations. It 
will also be interesting to find out using future XMM-Newton 
observations, whether the source makes a transition to a state 
which has a regular steep X-ray spectrum. Such an analysis 
will also shed light on whether or not the temporal property of 
the source, especially the observed QPO, depends on its un- 
usual spectrum. Search and analysis of other ULX that have 
similar hard spectra needs to be undertaken to verify whether 
M82 X-l is indeed a unique source. 
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